SMOOTHING EFFECT OF QUENCHED DISORDER 
ON POLYMER DEPINNING TRANSITIONS 



GIAMBATTISTA GIACOMIN AND FABIO LUCIO TONINELLI 

Abstract. We consider general disordered models of pinning of directed polymers on 
a defect line. This class contains in particular the (1 + l)-dimensional interface wetting 
model, the disordered Poland-Scheraga model of DNA denaturation and other (1 + d)- 
dimensional polymers in interaction with flat interfaces. We consider also the case of 
copolymers with adsorption at a selective interface. Under quite general conditions, 
these models are known to have a (de)localization transition at some critical line in the 
phase diagram. In this work we prove in particular that, as soon as disorder is present, 
the transition is at least of second order, in the sense that the free energy is differentiable 
at the critical line, so that the order parameter vanishes continuously at the transition. 
On the other hand, it is known that the corresponding non-disordered models can have 
a first order (de)localization transition, with a discontinuous first derivative. Our result 
shows therefore that the presence of the disorder has really a smoothing effect on the 
transition. The relation with the predictions based on the Harris criterion is discussed. 
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1. Introduction and models 

Quenched disorder is expected to smooth phase transitions in many situations. This 
is for instance the case of ferromagnetic spin systems in dimension d < 2 (if the spins 
are discrete) or in dimension d < 4 (if they have rotation symmetry): when a random 
magnetic field is present the Imry-Ma argument |2L)| . made rigorous by Aizenman and 
Wehr pP, implies that these models do not exhibit, at any temperature, the first order 
phase transition with associated spontaneous magnetization which characterizes the corre- 
sponding pure models. For the analogous phenomenon for SOS effective interface models 
see jS]. 

In the present work, under some conditions on the disorder distribution, we prove that 
a similar effect takes place in models of directed polymers in random media exhibiting a 
localization/delocalization transition on a defect line. Such a transition may be of first or 
higher order in the corresponding pure, i.e. non-disordered, cases and we show that, as 
soon as disorder is present, the transition is at least of second order. It is important to 
emphasize that the mechanism inducing the smoothing of the transition in our case is very 
different from the Imry-Ma one, and it is rather based on an estimate of the probability 
that the polymer visits rare but very favorable regions where the disorder produces a large 
positive fluctuation of the partition function. 
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We will consider mainly two classes of models: random pinning (or wetting) models and 
random copolymers at selective interfaces. In pinning models P| ll()| H!T |l24| l?7] the typical 
situation one has in mind is that of a directed path in (1 + d) dimensions, which receives a 
reward (or a penalty) at each intersection with a 1-dimensional region (a defect line, in the 
physical language) according to whether the charge present on the line at the intersection 
point is positive or negative. On the other hand the copolymer model, whose study was 
initiated in ^1] in the theoretical physical literature (see |23j and references therein for 
updated physics developments) and in in the mathematical one, aims at modeling a 

(1 + 1) -dimensional, directed heteropolymer containing both hydrophobic and hydrophilic 
components, in presence of an interface (the line S = in the notations of Section 11.2)1 
separating two solvents, situated in the upper and lower half-planes (S > and S < 0). 
One solvent favors the hydrophilic components and the other favors the hydrophobic ones, 
i.e., if the charge of the n-th monomer is positive (negative), this monomer tends to be in 
the upper (lower) half-plane. 

The main question one would like to answer is whether or not the interaction induces 
a localization of the polymer along the defect line or interface. 



1.1. Random pinning and wetting models. Let S := {S n } n=01 be a homogeneous 
process on an arbitrary set that contains a point 0, with So = and law P. For (3 > 0, h 6 
R and u> = {u> n }n=i,2,- - £ R N > one introduces the probability measure 



dP^_ {s) = exp ^ ^ _ h) lg \ (L1) 



dP V ' Z Na 



The choice of setting Sjy = is just for technical convenience, see Remark 1 1 . 1 1 b elow . Let 
us set To = and, for i G N := {1, 2, . . .}, n = inf {n > 7*j_i : S n = 0} if Tj_i < +co and 
Tj = +oo if T{-i = +oo. We assume that r := {ri} i is the sequence of partial sums of an 
IID sequence of random variables taking values in N U {oo} with discrete density K(-): 

K(n) = P(n = n). (1.2) 

This is of course the case if S is a Markov chain and for definiteness one should keep this 
case in mind. In order to avoid trivialities, we assume that K{oo) < 1 and for the model 
to be defined we assume also that there exists s € N such that P(SV = 0) > for every 
N G sN. Therefore, starting with p.ljl . we always implicitly assume that N £ sN, when 
iV is the length of the polymer. 

Moreover the cases we will consider are such that there exist a > 1 such that 

hminf^^>-a, (1.3) 

n— >oo log n 

and K{n) = if n sN. We advise the reader who feels uneasy with the weak requirement 
()1.3|) to focus on the case of K(n) behaving like n~ a , possibly times a slowly varying 
function (see Appendix [XJ. Note that, for a < 2, the return times to zero of S are not 
integrable. We stress that we have introduced s to account for the possible periodicity of S: 
of course the most natural example is that of the simple random walk on Z, P(5"j — = 
1) = p(Si - = -1) = 1/2 and {Si - S^i IID, for which (JEU) holds with s = 2 and 
a = 3/2 [111 Ch. III]. Another example is that of the simple random walk on Z d , d > 2: 
in this case, a = d/2. 
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We can rewrite Z^ jtJ , and in fact the model itself, in terms of the sequence r: the 
partition function in (jl.lj) is 

Z N ^ := E [exp (H N , U (S)) ; t Nn = N] , (1.4) 

with A/at := sup{i : Tj < N} and 

H N ,u( S ) = J2^~h). (1-5) 

i=l 

The disordered pinning model, which we consider here, is obtained assuming that the 
sequence oj is chosen as a typical realization of an IID sequence of random variables with 
law P, still denoted by u> = {u> n } n . We assume finiteness of exponential moments: 

E[exp(^i)] < oo, (1.6) 

for t € M and, without loss of generality, Eui = and EuJ = 1. Further assumptions on 
P will be formulated in Section where we state our main results. 

Under the above assumptions on the disorder the quenched free energy of the model 
exists, namely the limit 

F{/3,h) := lim ^]ogZ Nfi> , (1.7) 

N^co iv 

exists P( da;)-almost surely and in the L 1 (P) sense. The existence of this limit can be 
proven via standard super-additivity arguments based on Kingman's sub-additive ergodic 
theorem |22| (we refer for example to j2j for details). This approach yields also 
automatically the so called self-averaging property of the free energy, that is the fact that 
f(P, h) is non random. A simple but fundamental observation is that 

F(/U)>0. (1.8) 
The proof of such a result is elementary: 



lElog^o, > ^ElogE 



exp y^ifi^-h) );t 1 = N 



,i=l 



(1.9) 



= l(/3E[u; N ]-h) + ^logP(T 1 = N) 0, (1.10) 

where we have used assumption (jl.3|) . The proof of ()1.8|) suggests the following partition 
of the parameter space (or phase diagram): 

• The localized region: £ = h) : f(/3, h) > 0}; 

• The delocalized region: V = {(/3, h) : f(/3, h) = 0}. 

We set h c (P) := sup{/i : f(/3, h) > 0}, and we will call h c {j3) the critical point. Since 
f(/?, •) is not increasing and continuous, V = {((3, h) : h > h c ((3)}. It is rather easy to see 
that \h c (P)\ < oo, see e.g. |2j and [Tfij . 



Remark 1.1. It would be of course as natural to consider the model with partition 
function Z f N w := E [exp (Hn,u>{S))]- It is therefore worth to stress that by standard 
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arguments, see e.g. [TJI|, one sees that 



Zn,u> < Z f Nuj < max ^ " }^ exp (\(3uj n - h\) Z Nu; 
' " nesN: Kin) 
n<N v ' 

< — 777 r exp(|^a;Ar -h\)Z N(jJ , 

mm nesN: K(n) 

n<N 



(1.11) 



uniformly in N € sN and to. In particular, by ()1.3|) . the free energy is unaffected by the 
presence of the constraint Tj\f N = N: we stick to the constrained case because it simplifies 
some technical steps of the proofs. 



Remark 1.2. For ease of exposition we have used the generic denomination of pinning 
model, but our framework, as it is possibly clearer when the partition function is cast in 
the form (|1.4|) . includes a variety of models, like for example the Poland-Scheraga model 
of DNA denaturation, for which theoretical arguments on models with excluded volume 
interactions suggest a value of a larger than 2 |2J : note that whether the transition in the 
disordered case is of first or higher order is a crucial and controversial issue in the field, 
see for example [HI HIE]- We stress also that, since we can choose K(oo) > 0, it is rather 
easy to see that also the disordered (1 + 1) -dimensional interface wetting models jlUI I13| 
enter the general class we are considering. 

It is known that the nonrandom case (3 = is exactly solvable and that, according to 
the law of n, the transition at h c (0) = log(l — K(oo)) can be either of first or higher order. 
For completeness and to match our set-up, we give a quick self-contained analysis of this 
case in Appendix ^] 

1.2. Random copolymers at a selective interface. In the case of the copolymer 
model, the natural setting is to assume, in addition, that the state space of the process 
{S n } n is Z, that the law P is invariant under the transformation 5 — ► —S and that 
(Si — Si-i) £ {— 1,0, +1} for every i > 1. For instance, if {Si — Si_i}i is a sequence of 
IID variables then it is a classical result [El Ch. XII.7] that K(n) = c(l + o(l))n~ 3 / 2 , 
for large values of n £ sN (so a = 3/2): c is a constant that can be expressed in terms of 
P(S*i — So = 0) and of course s = 1 unless P(Si — So = 0) = 0. The copolymer model is 
defined introducing 

1 - 

n N ,US) = 75 5Z(/^n + h) sign(5 n ), (1.12) 

1 n=l 

with the convention that sign(0) = +1, and the corresponding partition function 

Z N := E [exp (Hn^S)) ; S N = 0] . (1.13) 

We may assume without loss of generality that both (3 and h are nonnegative. The factor 
1/2 in (|1.12|) is introduced just for convenience (see Section l3.2fl . As in Section ll.l[ the 
corresponding random model is obtained by choosing u as a realization of an IID sequence 
of centered random variables of unit variance and finite exponential moments. In analogy 
with (fT7|) and (fL8|) . the limit 

f(P,h) := lim ^\ogZ N ^ (1.14) 
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exists P( du;)-almost surely and in the L 1 (P) sense and one can prove as in (|1,9|) that 
f(P, h) := /(/?, h) — h/2 > 0. Therefore, also in this case we can partition the phase 
diagram into a localized and a delocalized phase, as 

• The localized region: £ = h) : f(/3, h) > 0}; 

• The delocalized region: V = {([5, h) : f(/3, h) = 0}. 

Also in this case, we set /i c (/?) := sup{/i : f(/3, /i) > 0} and we observe that V = {((3, h) : 



Many results have been proven for copolymers at selective interfaces, but they are 
almost always about the case of simple random walks. However they can be extended in a 
rather straightforward way to the general case we consider here (we omit the details also 
because they are not directly pertinent to the content of this paper). Above all we have 
that < h c {(3) < oo for every j3 > 0. Even more, the explicit bounds carry over to the 
general S we consider here: the upper bound is even independent of the choice of the 
law of S, while the lower bound [I] depends on a. 

We take this opportunity to stress also that various results about path behavior in the 
two regions are available, both for the pinning problem and for the copolymer, and, once 
again, mostly in the simple random walk case a = 3/2. For instance, it is known that in 
C long excursions of the walk from the line S = are exponentially suppressed and the 
fraction of sites where the polymer crosses the line remains nonzero in the thermodynamic 
limit (e.g., cf |16|I25| and references therein, and ^Ej)- On the other hand, in the interior of 
T> the number of intersections is, in a suitable sense, 0(log N) and, for the copolymer, 
the number of steps where sign(S' n ) = —1 is also O(logA^) [T7] . 



While the free energy can be proven to be infinitely differentiable with respect to all 
of its parameters in the region £, ^H], no results are available about its regularity at the 
critical point, apart from the obvious fact that F is continuous since it is convex. The 
result of the present paper partly fills this gap, showing that the transition is at least of 
second order, as soon as disorder is present. 

In order to state our main theorem, we need some further assumptions on the disorder 
variables to. We will consider two distinct cases: 

CI: Bounded random variables. The random variable lo\ is bounded, 



h > h c {(3)}. 



2. Smoothing of the depinning transition 



uji\ < M < oo. 




C2 



Unbounded continuous random variables. The law of oj\ has a density P(-) 
with respect to the Lebesgue measure on R, and there exists < R < oo 
such that 





in a neighborhood of x = 0. This is true in great generality whenever P(-) 
is positive, for example when the disorder is Gaussian and, more generally, 
whenever P(-) = exp(— V(-)), with V(-) a polynomial bounded below. 



Then, one has: 
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Theorem 2.1. Under condition CI or C2, both for the copolymer and for the pinning 
model, for every < j3 < oo there exists < c(/3) < oo ; possibly depending on P ; such that 
for every 1 < a < oo 

F(P,h) < ac(p)(h c (P) - hf (2.3) 

ifh<h c (p). 

Although the above result, coupled of course with (|1.8|) . seems to be in the same spirit 
as the rounding effect proven by Aizenman and Wehr U for the two-dimensional Random 
Field Ising Model, the physical mechanisms of smoothing are deeply different in the two 
cases. While is based on a rigorous version of the Imry-Ma argument |2L)| (i.e., a 
comparison between the effect of boundary conditions and of disorder fluctuations in the 
bulk due to the random magnetic field) in our case the boundary conditions play no role 
at all and everything is based on an energy-entropy argument inspired by [I]. 

Remark 2.2. It is important to observe that, as explained in Appendix El m the pin- 
ning case the deterministic model, j3 = 0, has a first order phase transition whenever 
*l2n<=N n -K( n ) < +°o 2 , in particular when a > 2. Theorem 12.11 therefore shows that 
the disorder has really a smoothing effect on the transition. But more than that is true: 
for a G (3/2,2), f(0, h) at h c (0) is strictly less regular than f(/3, /i) at h c {(3). In fact 
f(0,/i) > (/i c (0) - h) 2 - 5 for 5 e (0, (2a - 3)/(a - 1)) and small values of h c (0) - h > 
(for sharper results, see AppendixEJ). Notice that this is in agreement with the so-called 
Harris criterion |19j . which predicts that arbitrarily weak disorder modifies the nature of 
a second-order phase transition as soon as the critical exponent of the specific heat in the 
pure case is positive. In the present situation, this condition corresponds just to a > 3/2. 
The Harris criterion also predicts that the critical behavior does not change if a < 3/2, 
which is compatible with Theorem 12.11 Rigorous work connected to the Harris criterion, 
in the Ising model context, may be found in [7j. 

As one realizes easily from the proof of Theorem 12.11 given in Section |3 the constant 
c(/5) in (|2,3|) can be very large (of order 0((3~ 2 )) for f3 small. This is rather intuitive: 
for (5 — > one approaches the deterministic situation, where the transition can be of first 
order. 

Remark 2.3. In the theoretical physics literature, the (de) localization transition is claimed 
to be in some cases of order higher than two |28| . or even of infinite order |27| I23j. The 
method we present here, which is rather insensitive to the details of the model, does not 
allow to prove more than second order in general. It is likely that finer results require 
model-specific techniques. 

Remark 2.4. A generalized model: copolymers with adsorption. It is also possible to 
consider copolymer models with an additional pinning interaction |26j . as 

N N 

n N ,^(S) = (Jhun - hi) l Sn =o + 2 ]T(/^ n - h 2 ) sign(5 n ). (2.4) 

n=l n=l 

Here, both ui and oo are sequences of IID centered random variables with unit variance and 
finite exponential moments. In addition, one assumes u to be independent from u. This 
model corresponds to the situation where the interface between the two solvents is not 
neutral. While for simplicity we will not present details for this model, we sketch here what 
happens in this case. In analogy with the previous models, one can partition the phase 
diagram (i.e., the space of the parameters 0%, hi, 02, h^) into a localized and a delocalized 
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region, separated by a critical surface. In this case, Theorem 12.11 is easily generalized to 
give that the free energy f(/3i, hi, fh, h%) has continuous first derivatives with respect to 
/ii,/i2 when these parameters approach the critical surface from the localized region. 

3. Proof of the smoothing effect 

3.1. The pinning case. In this section we prove Theorem 12.11 for the pinning case, and 
in the next one we explain how the proof can be immediately extended to the copolymer 
model. 

The key idea, in analogy with |17j . is to introduce a new free energy where the fraction 
of sites where the polymer comes back to zero is fixed. In other words, recalling that 

M N = |{1 < n < N : S n = 0}| (3.1) 

one introduces, for m S [0, 1], 



(/>(/?, to) = lim lim — Elog Znu>(P; m, e) 

e\0 7V^oo N 



RSj&o ^ E1 ° gE {^ E "= lUlnlSn=01 (M N /N)e[ m -e, m+ ^S N = 



(3.2) 



Note that the limit is well defined, since Elog Zj^^(f3;m,e) is super-additive in N, thanks 
to the IID assumption on the increments of the sequence of return times r, and non- 
increasing for e \ 0. Therefore the limit 4> e (P, m) of (l/iV)E log Z]y^((3; to, e), as N — > oo, 
exists, as well as the second limit (j) e (P, m) \ 4>((3, m) as e \ 0. Notice moreover that (|3.2() 
holds also without taking the expectation, as for (1/-/V) log Zn^: in this case of course the 
limit N — > oo has to be taken in the P( do>)-a.s. sense. Moreover, it is immediate to realize 
that (p(f3, 0) = and that 

#,m)<F(A0)</3, (3.3) 

so that 4>{j3, m) is always bounded above. Finally, always thanks the IID property of the 
differences of successive return times to zero, it is easy to show that <fi({3, •) is concave: 

to) > x(p((3, mi) + (1 - x)(p((3, m 2 ) (3.4) 

if < x < 1 and m = xmi + (1 — x)m2- By exploiting the P(dw)~a.s. convergence of 
(1/N) log Zn jUJ ((3; to, e) to the nonrandom limit cf) e (f3,m) and the subsequent convergence 
for e \ 0, one deduces that f(/3, h) and 4>{(3, to) are related by a Legendre transform: 

F(j3,h)= sup (<j>(p,m)-hm). (3.5) 

me[0,l] 

In turn, this allows to identify h c ((3) in terms of (f>((3, •) as 

h c (P) = inf \ h : sup (<£(/?, m)-hm) = o\. (3.6) 

[ me[0,l] J 

The key technical step in the proof of Theorem 12.11 is the following: 

Theorem 3.1. Under condition CI or C2, for every < (3 < oo £/iere exists < C(/3) < 
oo such that, for every 1 < a < oo, 

009, to) - ft c C9)m < -^to 2 (3.7) 

a 

j/0 < TO < 1. 
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Proof of Theorem \2.1\ It is an immediate consequence of Theorem 13 . 1 1 and (|3.5|) . In fact, 
by (|3.7j) we have 

C(3) 

6(B, m)-hm< y -^m 2 + (h c (/3) - h)m, (3.8) 

a 

for every m € [0,1]. Taking the supremum over m on both sides of ()3.8|) . by (|3.5() we 
obtain 

F(P,h) < sup f-^lm 2 + (h c {0) - h)m) (3.9) 

m€[0,l] \ a J 

< sup f-^m 2 + (h c ((3) - h)m) = j^-AhM - h) 2 . 

Theorem 1^11 pinning 

We go now to the proof of Theorem 13.11 We will consider first the case in which lo\ 
satisfies condition C2, because it technically lighter. The two cases differ only in the 
first part of the proof (that is, up to Remark 13.21 below), where the probability of a rare 
event is estimated from below by changing the law of the disorder and by evaluating 
the corresponding relative entropy price. In the case of C2 it is sufficient to shift (i.e. 
to translate) the distribution of the disorder variables, and the relative entropy estimate 
implied by (|2,2j) fits well the rest of the proof. Under assumption CI, instead, we have to 
tilt the law of to and the arising expressions need to be re-worked, see Lemma 13.41 below, 
before stepping to the second part of the proof. 

Proof of Theorem 1,9. 1\ under assumption C2. Due to the concavity of 4>(f3, •), it is enough 
to prove (|3.7|) for < m < c\, with c\ = c\{j3) > 0, not depending on a. We define 

B 2 m 2 

7 (/3, m) := -<X/3, m) + h c {(3)m + c 2 - > 0, (3.10) 

a 

where C2 > is a constant depending only on P, which will be chosen later. We stress 
that the term containing C2 has been added simply because a priori one knows only that 
—<p(f3, m) + h c {(5)m > 0, cf. ()3.6|) . and it turns out to be technically practical to work with 
7(/3, m) > 0. For i G sN we define also 



wet': log Zt^(J3; m, e) - h c ((3)m£ > 7 (/3, m)^} . (3.11) 



Moreover for < m < c\ , we let P be the law obtained from P shifting the distribution of 
uji, • • • , uj£ so that 

,7(Am) 
(3m 

Note that 



E[uj] = 8 ■ (3.12) 



lim 7(/3, m)/m = 0, (3.13) 

since otherwise, by convexity of 7(/5, •), one has 7(/3, m) > em for for some e > and every 
m, that is f(/3, h c {(3) — e) = 0, which is in contrast with the definition of h c {(5). 

We now show that the event 

E := {u : (wi, . . . , u e ) E A i>m>£ } , (3.14) 
becomes P-typical for I large. 
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We first observe that, thanks to the constraint Ni/£ > m — e, and assuming that 
e < m/2, one has 

I (hog Z e Jfc m, e) - h c {/3)m) > 4 7 (/3, m) + E (hog Z e>w ((3; m, e) - h c {f3)m) .(3.15) 



J 

By ()3.2|) and (|3. lUf) . there exist £o(m) > and £o(e,m) € sN such that 

E (~ log Z l>ul (0; m, e) - h c {(3)r^j > 2 7 (/3, m), (3.16) 

for e < Eq, £ > £q. This in turn implies that 

F(E) >F (log Z et UP;m,e) - Elog Z^(/3; m, e) > - 7 (/3,m)/) , (3.17) 

which is greater than, say, 1/2 for £ sufficiently large, since 7 (/3, m) > and (recall 1)3.2)) 
and discussion following that formula) 

j\og Z^{l3-m,e) - ^Elog Z^{(5- m, e) ^ 0, (3.18) 

in P-probability. 

The price of shifting P to P is directly estimated by using the assumption C2, cf. 1)2.2(1 . 
and recalling 1)3.13)1 : assuming that m < ci(/3), with ci(/3) sufficiently small, one obtains 
the estimate 

H(F|P):= E feo g |W«(^\ (3.19) 



dP °dP|- \ /3m 
and by applying the relative entropy inequality 



we obtain 



for large 



log ( 5^ | > -J— (H(P|P) + e- 1 ) , (3.20) 

\p(£) y p(e) v / 

p:=P(E) > ^exp (-128^(^/3, m)/(/?m)) 2 /) , (3.21) 



Remark 3.2. Inequality 1)3.20)1 holds whenever the measures P,P are absolutely contin- 
uous with respect to each other and for every event E of nonzero measure. It is a simple 
consequence of Jensen inequality: since rlogr > — e _1 for every r > 0, one has 



um) =io gE (* 

\F(E) J V dP 



/ 


> E ^lo 


dP 


, dP 




log ——1 


dP 


& dP 



dP 

dip 



(3.22) 



1 / dP . dP \ 1 / / dP . dP 



E > E — log — + e 



-i 



F(E) \ dP d P / ~ $>(E) I \ dIP dP 
which is just (|3.20|) . 

We now apply an energy-entropy argument similar to that of |3j which, in the present 
case, roughly consists in selecting only those polymer trajectories which visit the rare 
stretches where the disorder configuration is such to produce a sufficiently large positive 
fluctuation of the partition function. Of course the precise definition of these rare stretches 
is directly related to the event E. This selection strategy gives a lower bound on the free 
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energy, which implies (|3.7|) , More precisely, we consider a system of length k£, with fcsN, 
£ G sN and we divide it into blocks Bj = {j£ + l,j£ + 2,--- , (j + !)£} of length £, with 
j = 0, • ■ • ,k — 1. For a given realization of u>, we denote by I(to) the ordered set of 
nonnegative integers 

I ( UJ ) = O'l, • • • . J|X(w)|} := {0 < J < — 1 : • • • ,tO(j+i)g) G ^, ra , £ }, (3.23) 

where it is understood that e < eo( m ) 5 ^ > ^o( m i £ ) an d m < ci as above, so that ()3.21j) 
holds. We bound the partition function Zn,w below by inserting in the average over the 
paths the constraint that Si ^ whenever i G Bj with j £ 1(oS), that Si = whenever 
i = j£ or i = (j + 1)£ with j G I(uj), and that, for every j G 

\{i G 5,- : 5, = 0}| r 

^ ^G[m-e,m + e]. (3.24) 

In this way, recalling the definition of At me and assuming that s < "f(f3,m)/(2\h c (f3)\), 




Figure 1. A typical trajectory contributing to the lower bound 113.251 . In this exam- 
ple k = 22, £ = 8, = {4, 10, 17}. Note that S„ / for n in a block Bj with 
j {4, 10, 17}, except at the boundary with a block Bj with j £ {4, 10, 17}, since by 
construction S n — at the boundaries of these blocks. Inside Bj, j £ , the path 
moves with relative freedom, but it is bound to touch the line S — approximately m£ 
times (see the text for the choice of m). 



one obtains 

llogZ^(/3,/ lc (/3))>i|ZH|^^ + llog K(L r ), (3.25) 

r=0,-,|Z(w)|: 
L r ^0 

where we recall that If (n) is the P-probability that first return to of an excursion of 
the free process occurs at step n, as in (|1.2|) . while L r 's, the (possibly vanishing) lengths 
of the excursions of the process between two blocks B,i,Bi> with i, i' G T(u>), are defined 
as 

L r = £\j r +i - jr ~ 1|, (3.26) 

with the convention that jo = —1, J|x(c^)|+l = see Figure^ Taking the expectation with 
respect to the disorder and using (|1.3j) . one obtains then 

-ElogZ k ^(f3,h c ((3)) > l^pp- -^E ^ l^ologi r; (3.27) 

r=0 
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for I sufficiently large. At this point, as in [I], one uses Jensen's inequality and the 
concavity of the logarithm to get 



iElog^ v (/?A(/3)) > ^f^P-^ 

> /VM ' > -2aE 



(|X(u/)| + l)log 



E 



|X(w)| T " 
r=0 n r 



k£ 



log 



ki 

\T(u)\ + l 



(3.28) 



Since the disorder variables in the distinct blocks {Bj}- are independent, the law of large 
numbers implies 



k 



p 



(3.29) 



dw)-a.s., so that, recalling (|3.21 
1 



= F(P,h c (p)) = lim -Elog Z H ,Mh c {/3)) 

k— »oo fix 

.7(W 2ap £ / 

> p — o r~ log ~ - p 

2 I p \ 



256aR 



(3m 



2a 



log(2£) 



(3.30) 



Since £ is arbitrary, one obtains 



j((3,m) > 



P 2 m 2 



(3.31) 

which is the desired inequality and l|3.7|) . provided that C2 in (|3.10f) satisfies C2 < (512i?) _1 . 

Theorc mpm C2 

Remark 3.3. It is easy to check that, in the Gaussian case uj\ ~ 7V(0, 1), Theorem 12.11 
holds with c(/3) = C3/3 -2 , C3 a suitable positive constant. Indeed, in this case the estimate 
(|2.2|) holds for every i£l (and R = 1/2) and therefore one can take c\ = 1 in the proof 
of Lemma EU so that (pOT)) implies (JEHI) with C(/3) = c 4 /3 2 . Of course, the same is true, 
up to constants, for every P such that inequality (|2.2I) holds uniformly for 

Proof of Theorem \S.l\ under assumption CI. The proof proceeds as in case C2, up to the 
definition of the law P, but in this case the law obtained by shifting P in general has an 
infinite entropy with respect to P. Therefore, in this case we define P rather by tilting the 
law of the first t variables: 



dP 1 ( v—-\ N 

— (u) = ? exp \uJ2"n 



n=l 



where u > will be chosen later and z = z(u) = Ee™ 1 . Let, for £ € sN, 



and observe that 



ij)t(u) = -E logZ £ J/3;m,s) 



^(0) = is log Z £ J (3; m,e). 



(3.32) 



(3.33) 



Then, one has 
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Lemma 3.4. There exist uq{(5) > and cq((3) > 0, possibly depending on F, such that for 
every < (3 < oo, 1 < a < oo the following holds: for every m £ (0, 1], if e < m/2 and 
< u < uq we have 

^i{u) - tp £ (0) > copmu. (3.34) 

Lemma 13.41 will be proven below. To proceed with the proof of Theorem 13. 11 we choose 
u = 4'y(P,m)/(Pmco) and notice that u is certainly smaller than no if m < c\{(3) with c\ 
sufficiently small (see (J3.13)) ). Then, choosing < m < ci and e < m/2, (j3.34j) implies 
that (|3.15|) is valid also in the present case. On the other hand, it is immediate to verify 
that (|3.19|) still holds, with R replaced by some R'((3, M) and m) by j((3, m)/co. The 
rest of the proof proceeds exactly as in the case C2. 

Theorem IO CI 

n— 

Proof of Lemma Note that 

1 1 



(3.35) 



i=l 



where ^ = £(u) = El<ji. The first term in the right-hand side of (|3.35j) can be rewritten 
(with some abuse of notation) as 



i=i 



where 



logZ e J(3; m,e) 



i=l 



i=l 



dy E^(l 5i=0 ) 
(ui -£) dy Ei Jl Si = ) 



E 



E^(-) 



' ex P ( Yh=X P^n ) 1 5 < =ol(A^/^)e[m-e,m+ £ ] 



(3.36) 



(3.37) 



and obviously the first term in (|3.36f) cancels the second one in the right-hand side of 
(|3.35|) . Next, observe that the following identity holds: 

Sy-y') e^(i Si=0 ) 

E^(ls»=o) 



1 + (eP{y-y') - l) E v (l Si=0 ) 



(3.38) 



Now recall f)2. lj) . so that it is sufficient to consider y and y' such that \y — y'\ < 2M, and 
from (|3,38|) we obtain 



-2PM 



< ^i,u;0-si=o) 



< e 



-2/3M 



^i=y 



(3.39) 



^i=y 



We can use this inequality to bound below the last term in (|3.36|) . We have in fact 



E 



dy E^(l 5i=0 ) 



>e- 2 ^r,EE^(l Si= 



>e-^ M nEB e>UJ (ls l=0 ) 



where n = n(u) = E tof, while 



E 



dyE^(l 5! =o) <Me 2 ^EE^(l 



s l =o) 



(3.40) 



^i=y 
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Therefore, recalling the constraint (Mi/£) £ [m — e, m + e] in the definition of P one 
has the following lower bound: 

duMu) > P(m - e)r?e- 4M/3 - (3(m + e)^Me 2(3M . (3.41) 

Now choose e < m/2 and notice that £ = u + 0(u 2 ) for u <C 1, while 77 = 1 + 0(u). 
Therefore, there exists uo((3,M) > such that, for < u < uq and for every m, the 
following holds: 

p -4M/3 e -4Af/3 

0«tM«) > Pm 2pmuMe 2 ^ M > f3m =: c Q f3m. (3.42) 

4 8 

An integration in u concludes the proof of (|3.34|) . 

Lemma 13. 41 



3.2. The copolymer case. In order to prove Theorem 12.11 for the copolymer model 
()1.12j) . it is convenient to start from the observation that one can rewrite the limit free 
energy as 



F(/3,h) = lim — ElogE 



exp 



N \ 

■J2(Pu n + h) A n 1^=0 

n=l / 



(3.43) 



where A r 



L sign(S'„)= 



Comparing this expression for the copolymer free energy with 



(jl.lj) . it is clear that in the present case the role of ls„=o is played by A n . The proof then 
proceeds exactly like in the pinning case, with the only differences that in the definition 
()3.2|) of (j>((3, m) the constraint (Mn/N) G [m — e, m + e] has to be replaced by 



|{1 < n < N : A n = 1}\ 



N 



£ [m — e, m + e] 



(3.44) 



and that, in the energy-entropy argument, the path is required to satisfy S% > (and not 
just Si ^ 0) whenever i € Bj with j £ This implies that K(L) in ()3.25j) has to be 

replaced by K(L)/2, which has the effect of adding a negative term of order 0(p/l) in the 
lower bound (J3.3U|) . which is negligible for I sufficiently large. 

Theorem 13^11 copolymer 



Appendix A. The non-disordered pinning model 

For (3 = 0, the free energy (|1.7j) may be identified explicitly with the following procedure. 
First we consider the equation 

^ K{n) exp(-bn) = exp(/i), (A.l) 

nGN 

and we look for a solution b > 0, which exists only if X^ngN-^( n ) > ex P(^)j m which case 
it is unique. Then if we set K{n) := exp(— h — bn)K(n), K{-) is a discrete probability 
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density and one can write 



N 



(A.2) 



E [exp (-hAf N ) ;t ^ n =N] = Y^ Yl II exp(-h)K(£i) 

J'=l (4,.../j)eN i :i=l 

AT j 

= exp (6iV) ]T ]I * =: ex P ( 6iV ) 

i=i (£i,...,^)gnj : »=i 

and one easily sees that Gat is the probability that the random walk which starts at 
and takes positive integer IID jumps with law K(-) hits the site N. It is a classical 
fundamental result of renewal theory that limjvG/v = l/SneN n -^( n ) d Ch. XIII]. 
This of course implies that b = f(0,/i). On the other hand, if (|A.1|) admits no positive 
solution, by proceeding as in ()A.1|) and by setting simply K(n) := exp(— h)K(n), so that 
K(-) is a sub-probability density, one easily sees that f(0, h) = 0. So equation equation 
(|A.1|) contains all the information about the free energy. 

Let us then observe that (|A.1|) has a positive solution if and only if h < log(l — K(oo)) 
and therefore h c (0) = log(l — K(oo)) < (h c (f3) being defined before Remark II, lj) . The 
behavior at criticality can be extracted from (jA.lj) in a rather straightforward way too, 
but of course we need to make precise the requirement on K{-) beyond the lower bound 
Q: 

• The case X := ^2 n ^nK{n) < oo. This is a necessary and sufficient condition for 
the transition to be of first order. More precisely, for £ E (0, oo] 

f (0, h) = eXp( ^ c(0)) (MO) -h) + o ((/i c (0) - h)) , for h / h c (0). (A.3) 
Formula (|A.3f) follows since by Dominated Convergence if S < oo then 

K ( n ) exp(-bn) = exp(/t c (0)) - 6S + o{b) (A.4) 

nGN 

for b \ 0, while by a direct estimate lini5\o X^neN K(n)[l — exp(— bn)] = oo if 
£ = oo. Note that the condition S < oo holds, in particular, in the case where 

K(n) = L{n)/n a (A.5) 

for large n E sN, with a > 2 and L(-) a function varying slowly at infinity, i.e., 
a positive function such that lim r ^oo L(xr) / L(r) = 1 for every x > (see j3j for 
more details on slowly varying functions). 

• The case £ = oo. We set K(n) = ^2J^ n+ i K(n) (by this we mean the sum over 
j € N such that j > n: K(oo) is not included in the sum) and assume that 

N 

J2K(j)=L(N)N 2 - a , 

3=1 

for some function L(-) which is slowly varying at infinity. This is true, in particular, 
in the case (TOl) . By the easy (Abelian) part of the classical Tauberian Theorem 
[El Ch.XIII.5, Theorem 2] we have that X^=o exp(-bn)K(n) = c'6 Q " 2 L(l/6)(l + 
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o(l)) as b \ 0, with d 



= d{a) > 0. Therefore 



CO 



OO 




(A.6) 



n=l 



n=0 



= K(0)-c'6 Q - 1 L(l/6)(l+o(l)), 



and 



F{0,h) = (h c (0) 



h) 



y^- 1 ) L(l/(h c (0)-h)), for h / h c (0) 



(A.7) 



with L(-) a slowly varying function (see (1.5.1) and Theorem 1.5.12]). It is 
therefore clear that the transition is of second order for a € (3/2, 2] (we emphasize, 
for the case a = 2, that we are assuming that £ = +oo) and it is of higher order for 
a < 3/2. The value a = 3/2 is borderline and the order of the transition depends 
then on the slowly varying function L(-). In the case of one dimensional symmetric 
random walks with IID increments taking values in {— 1,0, +1}, a = 3/2 and 
L(n) = c(l + o(l)) for n large, c a positive constant, and therefore the transition 
is really of second order and not higher. 



We would like to thank Bernard Derrida, Thomas Garel, Cecile Monthus and David 
Mukamel for interesting discussions. This research has been conducted in the framework 
of the GIP-ANR project JC05_42461 (POLINTBIO). 



[1] M. Aizenman and J. Wehr, Rounding effects of quenched randomness on first-order phase transitions, 

Comm. Math. Phys. 130 (1990), 489-528. 
[2] K. S. Alexander and V. Sidoravicius, Pinning of polymers and interfaces by random potentials, preprint 

(2005). Available on: arXiv.org e-Print archive: math.PR/0501028 
[3] N. H. Bingham, C. M. Goldie and J. L. Teugels, Regular Variation, Cambridge University Press, 

Cambridge, 1987. 

[4] T. Bodineau and G. Giacomin, On the localization transition of random copolymers near selective 

interfaces, J. Statist. Phys. 117 (2004), 801-818. 
[5] E. Bolthausen and F. den Hollander, Localization transition for a polymer near an interface, Ann. 

Probab. 25 (1997), 1334-1366. 
[6] A. Bovier and C. Kiilske, There are no nice interfaces in (2 + 1) -dimensional SOS models in random 

media, J. Statist. Phys. 83 (1996), 751-759. 
[7] J. T. Chayes, L. Chayes, D. S. Fisher and T. Spencer, Correlation Length Bounds for Disordered Ising 

Ferromagnets, Commun. Math. Phys. 120, 501-523 (1989). 
[8] B. Coluzzi, Numerical study on a disordered model for DNA denaturation transition, 

|arXiv:cond^m at/0504080 /l (2005). 
[9] D. Cule and T. Hwa, Denaturation of heterogeneous DNA, Phys. Rev. Lett. 79 (1997), 2375-2378. 
[10] B. Derrida, V. Hakim and J. Vannimenius, Effect of disorder on two-dimensional wetting, J. Statist. 

Phys. 66 (1992), 1189-1213. 
[11] W. Feller, An introduction to probability theory and its applications, Vol. I, Third edition, John Wiley 

& Sons, Inc., New York-London-Sydney, 1968. 
[12] W. Feller, An introduction to probability theory and its applications, Vol. II, Second edition, John 

Wiley & Sons, Inc., New York-London-Sydney, 1971. 
[13] G. Forgacs, J. M. Luck, Th. M. Nieuwenhuizen and H. Orland, Wetting of a Disordered Substrate: 

Exact Critical behavior in Two Dimensions, Phys. Rev. Lett. 57 (1986), 2184-2187. 
[14] T. Garel, D. A. Huse, S. Leibler and H. Orland, Localization transition of random chains at interfaces, 
Europhys. Lett. 8 (1989), 9-13. 



ACKNOWLEDGMENTS 



References 



16 



GIAMBATTISTA GIACOMIN AND FABIO LUCIO TONINELLI 



[15] T. Garel and C. Monthus, Numerical study of the disordered Poland-Scheraga model of DNA denat- 

uration, J. Stat. Mech., Theory and Experiments (2005), P06004. 
[16] G. Giacomin, Localization phenomena in random polymer models, preprint (2004). 

Available online: http:/ /www. proba.jussieu.fr/pageperso/giacomin/pub/publicat. html 
[17] G. Giacomin and F. L. Toninelli, Estimates on path derealization for copolymers at selective interfaces, 

Probab. Theor. Rel. Fields 133 (2005), 464-482. 
[18] G. Giacomin and F. L. Toninelli, The localized phase of disordered copolymers with adsorption, preprint 

(2005). Available on: arXiv.org e-Print archive: math.PR/0510047 
[19] A. B. Harris, Effect of random defects on the critical behaviour of Ising models, J. Phys. C 7 (1974), 

1671-1692. 

[20] Y. Imry and S.-K. Ma, Random-Field Instability of the Ordered State of Continuous Symmetry, Phys. 

Rev. Lett. 35 (1975), 1399-1401. 
[21] Y. Kafri, D. Mukamel and L. Peliti, Why is the DNA denaturation transition first order?, Phys. Rev. 

Lett. 85 (2000), 4988-4991. 
[22] J. F. C. Kingman, Subadditive ergodic theory, Ann. Probab. 1 (1973), 882-909. 

[23] C. Monthus, On the localization of random heteropolymers at the interface between two selective sol- 
vents, Eur. Phys. J. B 13 (2000), 111-130. 

[24] N. Petrelis, Polymer pinning at an interface, preprint (2005). Available on: arXiv.org e-Print archive: 
|math.PR/0504464| 

[25] Ya. G. Sinai, A random walk with a random potential, Theory Probab. Appl. 38 (1993), 382-385. 
[26] C. E. Soteros and S. G. Whittington, The statistical mechanics of random copolymers, J. Phys. A: 

Math. Gen. 37 (2004), R279-R325. 
[27] L.-H. Tang and H. Chate, Rare-Event Induced Binding Transition of Heteropolymers, Phys. Rev. 

Lett. 86 (2001), 830-833. 

[28] A. Trovato and A. Maritan, A variational approach to the localization transition of heteropolymers at 
interfaces, Europhys. Lett. 46 (1999), 301-306. 

Laboratoire de Probabilities de P 6 & 7 (CNRS U.M.R. 7599) and Universite Paris 7 - Denis 
Diderot, U.F.R. Mathematiques, Case 7012, 2 place Jussieu, 75251 Paris cedex 05, France 
Home page: http: //www .proba. jussieu. fr/pageperso/giacomin/GBpage .html 
E-mail address: giacomin@math.jussieu.fr 

Laboratoire de Physique, ENS Lyon (CNRS U.M.R. 5672), 46 Allee d'Italie, 69364 Lyon 
cedex 07, France 

Home page: |http : //perso .e ns-lyon . f r/f abio -Tucio . toninelli | 
E-mail address: fltonine@ens-lyon.fr 



